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the Need for AGN Heating 
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Abstract. We present hydrodynamical A^-body simulations of clusters of galaxies with feedback 
taken from semi-analytic models of galaxy formation. The advantage of this technique is that the 
source of feedback in our simulations is a population of galaxies that closely resembles that found 
in the real universe. We demonstrate that, to achieve the high entropy levels found in clusters, active 
galactic nuclei must inject a large fraction of their energy into the intergalactic/intracluster media 
| throughout the growth period of the central black hole. These simulations reinforce the argument of 

Bower et al. [1], who arrived at the same conclusion on the basis of purely semi-analytic reasoning. 
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INTRODUCTION 



X-ray observations of groups and clusters of galaxies allow us to probe the physical 



properties of the hot, diffuse intracluster medium (ICM). In the self-similar collapse 
scenario, where the ICM is heated solely by gravitational processes, we expect the X- 
ray luminosity Lx of clusters to scale with gas temperature T as Lx °= T 2 . However, the 
observed Lx-T relation is much steeper than predicted, with Lx °= r 2 - 5 ~ 3 at T > 2 keV, 
becoming steeper still at group scales T < 2 keV. This deficit in X-ray luminosity is due 
to an excess of entropy 1 in cluster cores. The source of this excess entropy is likely to 
be a combination of non-gravitational cooling and heating processes. 

The most obvious sources of non-gravitational heating are Type II supernovae (SNe) 
and active galactic nuclei (AGN). Cutting-edge hydrodynamical Af-body simulations of 
galaxy groups and clusters attempt to couple cooling, star formation and black hole 
(BH) growth with associated feedback from SNe and AGN (e.g. [2]). Although in 
their infancy, such simulations have already yielded encouraging results. For example, 
Puchwein et al. [3] resimulated a sample of 21 groups and clusters with the model of 
Sijacki et al. [2], obtaining a mean Lx-T relation and halo gas fractions that are in good 
agreement with observational data over a wide range of mass scales. However, the stellar 
fraction within the virial radii, r v - ir , of their objects appears larger than observations 
suggest, implying that too much gas has cooled, even with stellar and AGN feedback. 



We define entropy as S = kT jn\ , where k is Boltzmann's constant, n e is the electron number density 
and 7 = 5/3 is the ratio of specific heats for a monoatomic ideal gas. 



In this work we pursue a different, but complementary, approach to the theoretical 
study of galaxy groups and clusters. Instead of undertaking fully self-consistent hydro- 
dynamical simulations, we investigate what current semi-analytic models (SAMs) of 
galaxy formation predict for the thermodynamical properties of the ICM. These mod- 
els have already proved capable of explaining many key observational properties of real 
galaxies. Our goal is to extend the predictive power of SAMs to group and cluster scales 
by coupling them to cosmological hydrodynamical /V-body simulations. The main ben- 
efit of this hybrid approach is that feedback in our simulations is guaranteed to originate 
from a realistic galaxy population. This is generally not the case in self-consistent hy- 
drodynamical simulations. 

THE NUMERICAL MODEL 

Or hybrid technique consists of three components which we now describe in turn. For 
full details of the modelling process, we refer the reader to Short and Thomas [4]. 

A dark matter simulation. We adopt a spatially-flat ACDM cosmological model with 
the same parameters as the Millennium simulation [5]. Initial conditions (ICs) were 
generated for a cubic volume of side length L = 125/i 1 Mpc containing A^dm = 540 3 
dark matter (DM) particles of mass moM = 8.61 x 10 8 /i 1 M . We evolved these ICs to 
z = using the Af-body/SPH code GADGET-2, storing particle data at the 64 output 
redshifts of the Millennium simulation. DM halos are identified on the fly with the 
friends-of-friends (FOF) algorithm and bound substructures orbiting within FOF halos 
are found in post-processing with the SUBFIND algorithm. Halo merger trees are then 
constructed as in [5]. 

Feedback from a semi-analytic model of galaxy formation. A galaxy catalogue is 
generated for our DM simulation by applying the L-Galaxies SAM described by Croton 
et al. [6] and De Lucia and Blaizot [7] to the merger trees. For each galaxy in the 
catalogue, we use its merger tree to compute the change in stellar mass, AM*, and 
mass accreted by the central BH, AMbh, between successive model outputs. The energy 
imparted to the ICM by Type II SNe is then calculated from AM* via Eq. (20) of Croton 
et al. [6]. To determine the amount of heat energy input into the ICM by AGN, we do not 
use the L-Galaxies AGN feedback prescription, but instead use the scheme employed in 
the version of the SAM GALFORM presented by Bower et al. [1]. The reason for this 
choice is discussed in Short and Thomas [4]. The energy transferred to intracluster gas 
by AGN then follows from AMbh by using Eq. (3) of Bower et al. [1]. 

A hydrodynamical simulation. The ICs for our DM simulation are modified by 
adding gas particles with zero gravitational mass. This ensures the DM distribution will 
be unaffected by the baryons, so the halo merger trees will be the same. A modified 
version of GADGET-2 is used to evolve these ICs. Whenever an SPH calculation is 
required, we assign the gas particles their true mass (assuming Qb,o = 0.045), so that 
gas properties are computed correctly. Once an output redshift is reached, temporary 
'galaxy' particles are introduced throughout the simulation volume at positions specified 



by the SAM galaxy catalogue. For each galaxy particle, we know the change in stellar 
mass and energy released by SNe/AGN since the last output (see above). We use this 
information to form stars and heat gas in the vicinity of each galaxy as detailed in Short 
and Thomas [4]. The galaxy particles are then removed and the simulation proceeds 
until the next output time, when the process is repeated. Note that we choose to neglect 
cooling processes in our simulation (see below). 

RESULTS AND DISCUSSION 

Two key observables of galaxy groups and clusters are the Lx-T relation and the halo 
gas fraction. We now investigate whether our model can explain the observational data. 

The X-ray luminosity-temperature relation. The Lx-T relation obtained from our 
hybrid simulation at z = is shown in the left panel of Fig. 1. It is clear that the slope 
and normalisation of our relation is consistent with observations at all mass scales. In 
particular, we see the same steepening on group scales. This is because AGN heating is 
more efficient at driving X-ray emitting gas from the central regions of low-mass halos, 
reducing the gas density and thus X-ray luminosity. Our results also exhibit a variation 
in scatter about the mean relation that is similar to the data for T < 3 keV. Scatter about 
the low-temperature end of the Lx-T relation is due to the variety of merger histories 
of groups. Note that there seems to be more scatter towards the lower-luminosity edge 
of the observed relation than the upper-luminosity edge. The reason for this is that we 
cannot produce systems with a highly X-ray luminous cool core since gas cooling is not 
incorporated in our simulation. 

Halo gas fractions. The halo gas fractions / gas of our groups and clusters are plotted 
as a function of emission- weighted temperature in the right panel of Fig. 1. It is clear 
that there is a broad agreement with the observational data, with our results exhibiting 
a comparable amount of scatter. As in the data, we see a rapid decline in gas fraction at 
lower temperatures. This is because AGN feedback is effective at removing gas from the 
central regions of groups since they have a much shallower potential well than massive 
clusters. In addition, the stellar fraction within r wn is, on average, approximately 9% for 
our objects, agreeing with observations (e.g. [8, 9]). This is to be expected since star 
formation in our simulation is driven by a SAM which has been tuned to reproduce the 
cosmic star formation history. By contrast, self-consistent hydrodynamical simulations 
typically predict 20 — 50% of the baryons within r v i r are locked-up in stars (e.g. [3, 8]). 

SUMMARY AND CONCLUSION 

In this work we set out to extend the predictive power of current SAMs of galaxy 
formation by investigating the effect of energy feedback from model galaxies on the 
properties of intracluster gas. To achieve this objective we have employed a novel hybrid 
technique in which a SAM is coupled to a non-radiative hydrodynamical simulation, thus 
guaranteeing that the source of feedback is a realistic galaxy population. 
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FIGURE 1. Bolometric X-ray luminosity (left panel) and halo gas fraction (right panel) as a function 
of emission-weighted temperature. All X-ray properties are computed within rsoo- Our simulated groups 
and clusters are shown by open black circles and several datasets from X-ray observations are also shown 
for comparison. In the left panel, the upper and lower solid black lines are best fit relations obtained from 
non-radiative and preheating simulations, respectively. In the right panel, the horizontal dashed line is the 
mean cosmic baryon fraction for the cosmological model we have adopted. 



Our main conclusion is that a large energy input from AGN (on average, 35% of 
the available rest mass energy O.IMbhc 2 ) is required over the entire formation history 
of halos in order to reproduce the observed Lx-T relation and halo gas fractions. This 
supports the conclusion of Bower et al. [1] derived using purely semi-analytic reasoning. 
Indeed, the relation we have obtained agrees well with observations on all mass scales 
probed, displaying a similar degree of scatter. This is also true of the halo gas and stellar 
fractions of our simulated groups and clusters. 

In future work we plan to self-consistently incorporate cooling into our hybrid ap- 
proach. We note that the inclusion of gas cooling could only lead to a decrease in the 
entropy of intracluster gas, which would reinforce our conclusion with regard to the high 
degree of AGN feedback that is required to explain the excess entropy found in clusters. 
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